Neutron Scattering and the B\ g Phonon in the Cuprates 
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The momentum dependent lineshape of the out-of-phase oxygen vibration as measured in recent neu- 
tron scattering measurements is investigated. Starting from a microscopic coupling of the phonon 
vibration to a local crystal field, the phonon lineshift and broadening is calculated as a function of 
transfered momentum in the superconducting state of YBa2Cu307. For the first time, the absolute 
magnitude of the shift is obtained as a function of momentum and is found to be in excellent agree- 
ment with experiment. It is shown that the anisotropy of the density of states, superconducting 
energy gap, and the electron-phonon coupling are all crucial in order to explain these experiments. 
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Recently attention has focused on the observance of a 
resonance at an energy near 41 meV in the neutron scat- 
tering cross section obtained in both optimally doped and 
underdoped YBa 2 Cu307_5 |^,^]. It has been speculated 
that this resonance feature might hold the key to the pair- 
ing mechanism which underlies superconductivity in the 
cuprate materials. As such, this feature has been lavished 
with attention by many various theoretical proposals || . 
Earlier measurements had suggested that this resonance 
was related to a particular phonon vibration which oc- 
curs at nearly the same energy. This phonon consists 
of the out of phase c— axis vibration of the two oxygen 
0(2) and 0(3) atoms in the Cu02 plane. However it was 
later clarified via kinematic analysis and polarized neu- 
tron studies that a magnetic and phononic contribution 
are both present and could be separated. Both peaks 
could be tracked as a function of temperature as well as 
scattered neutron momentum. The phonon peak (2| has 
received a much smaller amount of attention in compar- 
ison with the magnetic resonance. 

Earlier Raman || and neutron scattering J5| experi- 
ments have studied this particular phonon vibration in 
detail. Of all the Raman-active phonons measured in 
tetragonal superconductors, only this phonon does not 
transform according to the full symmetry (Ai g ) of the 
lattice. Since this phonon obeys its own selection rules 
(Big) it can be unambiguously identified. From its Fano 
profile, the magnitude of the electron-phonon (e-ph) cou- 
pling seems to be particularly large for this phonon com- 
pared to other phonons measured in Raman scattering. 
Indeed from the Raman point of view, this phonon has 
also been lavished with attention |] ^ . In Ref. ||] , a phe- 
nomenological form for the coupling was assumed and 
thus predictions concerning the absolute magnitude of 
the lineshift could not be obtained. Moreover, a simpli- 
fied momentum dependence of the coupling was assumed. 
Only Refs. @,| and @ specified a microscopic source of 
the coupling. In Refs. |7]|| the coupling was a result of 
the buckling of the Cu-0 plane. On the other hand, Ref. 



H has addressed this coupling as arising through a local 
crystal electric field E which breaks the inversion sym- 
metry locally around the Cu-0 plane |10|. Here the ab- 
solute value of the coupling estimated from calculations 
of the electric field gave good agreement to the results 
obtained from Raman scattering studies. Neither theo- 
ries made predictions about the momentum dependence 
of the phonon lineshape. 

Ref. Q showed that in the normal state the B\ g 
phonon does not show any appreciable dependence on the 
in-plane neutron's scattered momentum q. However in 
the superconducting state the phonon softened consider- 
ably for q = 0, in agreement with Raman measurements, 
while the effect became rapidly less pronounced with in- 
creasing q. Along the (1, 1,0) direction in the Brillouin 
Zone (BZ), the phonon hardened faster with increasing 
q than along the (1,0,0) direction. While these findings 
were related to anisotropics of the energy gap, density 
of states, and the e-ph coupling, an unambiguous de- 
termination of the importance of each effect could not 
be made. We therefore refocus attention on the phonon 
contribution to neutron scattering in an effort to provide 
a theoretical framework to better understand the contri- 
butions from various anisotropics. Using a microscopic 
coupling theory developed in the momentum depen- 
dence of the phonon lineshape will be used to determine 
the interplay of the anisotropies of the energy gap, den- 
sity of states, and the e-ph coupling. For the first time, 
an estimate for the magnitude of the lineshape changes 
as a function of momentum can be obtained. 

Our starting point is a three-band model for the Cu02 
plane with Cu-0 hopping amplitude t and 0-0 hopping 
amplitude t': 

Ho = sJ2 b i« b ^ + 1 { p 8 b l° a n,6,* + h - c -] 

n >< 7 (6,6') 
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where b^ n a creates an electron with spin a at a cop- 
per lattice site n, while a n g annihilates an electron at 
one of the neighboring oxygen sites n + 8/2 determined 
by the unit vector 8 assuming the four values, (±1,0) 
and (0, ±1). An oxygen atom between the two copper 
atoms at sites n and n + 8 is labeled by either (n, 8) or 
(n + 8,-8). As in |J, e = Ed — E p is the difference 
of the Cu and O site energies and Pg = ±1 depending 
on whether the orbitals (with real wavefunctions) have 
the same or opposite sign at the overlap region. Assum- 
ing Cu d x 2_ y 2 and O p orbitals P _g = —Pg, and we 
can choose P(i,o) = 1 an d P(o.i) = ~ 1- Lastly, Pj ^, 
denotes the overlap sign between an oxygen orbital at 
site n + 8/2 with a neighboring oxygen orbital at site 
n + 8/2. The 0-0 hopping is needed to give the right 
curvature and centering of the observed Fermi surface. 
By our above convention these overlaps take the values 

K.y = p -x,- y = i'-fx-y = F -x.y = -I) respectively. 
After Fourier transforming, the Hamiltonian now reads 



Ek,,C. where 



" (2) 
+{ ib t,A a x,k,<7 t xO < -) - a y>k ,aty{k)] + h.c.} 
+t'(k)[<4 k . CT ay,k,o- + h.c], 



with the prefactors t a (k) = 2t sin(afc Q /2), t'(k) = 
—4:t'sm(ak x /2)s'm(ak y /2). We can then diagonalize Eq. 

@ : H l,a = E/3 e /3( k )4,k, CT d /3.k,a, where (3 is +,-, and 
for the antibonding, bonding, and nonbinding bands, 
respectively. We consider only a reduced one band model 
appropriate near half filling and take only the upper band 
into account |ll[] . At this point the same procedure used 
in Ref. |)| can be carried through. The electric fields at 
the 0(2) and 0(3) sites couple linearly to their c-axis ion 
displacements leading to the following e-ph coupling 



Uri- 
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eEz ■ ^{113,(011)04 XCT a r 

n,<T 



u v (on)oi ) y j(r o I1 ,y, .}, 



(3) 



where E is the c— axis aligned electric field at the 0(2) 
and 0(3) sites, neglecting orthorhombic distortions. This 
can be rewritten in a momentum representation as 



C(k)^(k - q)e-^°/ 2 (l + e-** ) 
-0*(k)c^(k - q)e-^°/ 2 (l + e- l ^ a ) 



Here Mq is the oxygen mass and lo Bi ~ 41 meV is 
the phonon frequency. The <fi functions result from the 
diagonalization of the three band model and are given as 

<t>*M = 5-[*»,tf(k) - i'(k)^(k)/e(k)] (6) 

with the normalization factor 
1 



A k = 



6(k) 



[( £ (k) 2 -t'(k) 2 ) 2 + (6(k)^(k) 



t y (k)t'(k)) 2 + (e(k)t y (k)~t x (k)t'(k)) 



2ll/2 



(7) 



Lastly the exponential factors and the other normaliza- 
tion N(q) comes from the phonon displacements. The 
eigenvalue and momentum dependent eigenvector for the 
Big ionic vibration is solved for the simple harmonic 
model of 3 atom unit cell of Cu coupled via identical 
spring constants K to both the 0(2) and 0(3) atoms 
under a uniform tension. Besides an acoustic mode 
and an in-phase 0(2)-0(3) optical mode (Ai g ), the B lg 
mode is obtained with the dispersionless energy eigen- 
value uj Bi = 2K/Mq, and whose eigenvector is 



€^' J, (q) = ±2(1 + e-^ a )/^N^>, 
N(q) = 4[cos 2 (q x a/2) + cos 2 (q y a/2)}. 



(8) 



The eigenvectors from this simplified model give an ad- 
equate agreement to the eigenvectors obtained via more 
complex lattice dynamical calculations . 

The e-ph coupling renormalizes the bare phonon prop- 
agator Dg (q, a;) = 2u Blg /[uj 2 — io 1 Bi ] such that the 
retarded propagator becomes 



D Blg (q, co + iS) = 



l+D B Jq,Lu)n Blg (q,Lj + tS)' 
where Tl Blg is the q and u> dependent polarization 



n Blg (q,o; + ^) = ^| 5 (k,q) |' 



(9) 
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where c q annihilates the B\ g phonon mode with wave 
vector q and g(k, q) is the coupling constant of the mode 
to an electron of wave vector k. 

The coupling constant for the q = B± g phonon was 
evaluated in ||. These results can be generalized for 
finite q as 



;(k,q) = eE z 



2M N{q)uj Blg 



(5) 



x^a + (k,q)[/(£; k )-/(£ k+q )] 

1 1 



UJ + l5 + £k - i?k+q CO + iS - E k + £"k+q 

a_(k,q)[l-/(£ k )-/(£ k+q )] 

1 1 



oj + iS — E\r — E] 



k+q 



w + iS + Eu + E } 



k+q 



where E^ = yCjf-i- A 2 (k), £ k = e(k) — [i is the energy 
dispersion measured from the chemical potential /i, / is 
the Fermi function, and a±(k, q) are coherence factors 
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FIG. 1. Real and Imaginary parts of the B\ g phonon self 
energy as a function of momentum transfer q along the BZ 
diagonal for T = T c /2. The solid lines in the figures corre- 
spond to the d x 2_ y 2 energy gap and the long dashed lines 
correspond to the | d x 2_ y 2 | energy gap. The left panels rep- 
resent the evaluation of Eqs. jlC| , |ll| ) using the expression for 
the coupling constant g, Eq. (0)7 while the right panels corre- 
spond to using a momentum-independent coupling constant 
equal to t. 



«±(k,q) = 1± 



eiA+ q -A(k)A(k + q) 



(11) 



This expression is identical to the one used in Ref. 
||. However, all other uses of Eqs. (|l^,^l|) have used 
a momentum independent coupling constant and have 
not specified a microscopic mechanism for the electron- 
phonon interaction. We remark here this misses crucial 
information of the interplay of symmetry of the phonon 
vibration and its coupling to the electronic system. First 
of all, the use of k- independent coupling constant in 
Eqs. ( |i^ , |TTl ) corresponds to a particular case of a fully 
symmetric phonon (A\ g ), transforming as the identity 
representation. As a consequence of the long wavelength 
Coulomb interaction, the q = self energy would be 
completely screened, leading to no net linewidth change 
for q = phonons, as measured e.g. in Raman scat- 
tering measurements P,|l3[ . If a fully symmetric phonon 
has a momentum dependence, it still will be at least par- 
tially screened by the long wavelength Coulomb interac- 
tion. However, there is no screening for phonons of other 
symmetry. While it appears that certain A\ g phonons 
in the cuprates are weakly temperature dependent and 
may indeed have their coupling to the electron system 
screened out, the Ba A\ g phonon as well as the O B lg 
vibration show distinctive shifts as a function of temper- 
ature. This can only come about if the phonons have a 
non-trivial momentum dependent coupling constant and 
therefore the use of a momentum independent coupling 
constant is inappropriate for these phonons Q . Secondly, 
the neglect of the momentum dependence of the coupling 
constant overestimates the size of the B\ g phonon broad- 
ening shift for q = 0. As with the electronic contribu- 
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FIG. 2. Real and Imaginary parts of the Bi g phonon self 
energy as a function of momentum transfer q along the BZ 
axis for T — T c /2. Both gap choices yield identical results. 
The solid lines in the figures correspond to evaluating Eqs. 
ll| ) using the expression for the coupling constant g, Eq. 
while the dashed lines correspond to using a momen- 
tum-independent coupling constant equal to t. 



tion to the Raman response |1J] , the coupling of the gap 
and the vertex which have the same symmetry leads to 
a reduction of the phonon broadening in the supercon- 
ducting state ~ (w/2A) 3 compared to the case when the 
gap and vertex are of different symmetries ~ uj/A [ pi| |. 
Lastly, the momentum dependence is not separable into 
two functions of k and q respectively, as implied in Ref. 

Here we consider either an anisotropic s— wave or 
d x 2_ y 2 energy gap: A s (k) = Ao, | [cos (Ala; a) — cos(fcj,a)] | 
/2, Arf(k) = Ao[cos(k x a) — cos(k y a)]/2. While we be- 
lieve that there is strong evidence that the energy gap 
in the cuprates has d x 2_ y 2 symmetry, we consider the 
s— wave case for illustrative purposes showing the effect 
of a sign change of the energy gap. In the following 
we choose e — leV, t = 1.6eV, t'/t = 0.45 (which are 
similar to the values chosen in Ref. ||), Ao = 30meV 
and 2Ao/fcsX' c = 8 for both energy gap choices, and 
oj = 41meV. We adjust the chemical potential to yield 
a filling (n) = 0.875 for both spins, and set T = T c /2. 

In Fig. (1) we plot the real and imaginary parts of 
n as a function of transfered momentum q along the 
(1,1,0) direction for a d x 2_ y 2 and an anisotropic s— wave 
superconductor, while Fig. (2) plots these same functions 
for q along the (1,0,0) direction. We have normalized 
both the real and imaginary parts of n to their q = 
values. The three factors governing the q dependence 
of the B\ g phonon's lineshape are: the joint density of 
states for phonon scattering, the coherence factors, and 
the e-ph coupling. 

We first consider features of the the polarization n 
which are only due to kinematic constraints and band 
structure. Starting from q — 0, where Cooper pairs are 
broken equally around the Fermi surface (FS), the pair- 
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FIG. 3. The momentum dependence of the coupling con- 
stant | g(k, q = 0) | 2 as a function of momentum in the BZ. 




FIG. 4. The momentum dependence of the coupling con- 
stant | <?(k,q = w/a,0) 2 as a function of momentum in the 
BZ. 



breaking becomes less resonant as q increases until two 
regions of the FS can be connected. At the same time the 
scattering processes at low temperatures which are most 
prominent are those which can most closely satisfy the 
condition that uj =1 A(k) | + | A(k + q) | coming from 
last term in Eq. (|l0|). Both of these restrictions become 
less and less obeyed with increasing q and lead to a very 
rapid drop off of II away from q = in either direction as 
seen in both figures regardless of the energy gap symme- 
try. Moreover, as q approaches the zone boundary along 
the diagonal (Fig. 1), both the real and imaginary parts 
of II rise since regions of the FS with large density of 
states (flat bands) are connected by momentum transfer 
of (it /a, ir /a). 

Now if we consider the anisotropy of the energy gap, 
we note that for the anisotropic s— wave case this rise 
is more pronounced than the d— wave case due to the 
change of sign of the gap in the coherence factor a- for 
the d— wave gap. The rapid fall off of the the data in 
Ref. Q gives support for d x 2_ y 2 pairing. As we would 
expect, there is no difference between II evaluated for 
the two gaps for q transfer along the (1,0,0) direction 
since the coherence factor is the same for both energy 
gaps for scattering near the FS. 

We now consider the momentum dependence of the 
coupling constant. In the Fig. 1 (left panel) and Fig. 
2 the solid line represents II calculated with the cou- 
pling constant included while the right panel of Fig. 1 
and the dashed lines in Fig. 2 are obtained by neglect- 
ing the k, q dependence of the coupling constant. We 
see that the results for both q directions are strongly af- 
fected by the coupling constant and its inclusion cannot 
be neglected. In all cases the coupling constant leads to 
a more rapid fall off of the polarization with increasing q. 
This is due to the strong momentum dependence of the 
coupling constant <?(k, q). For q = the vertex varies as 
~ [cos^a) — cos(fc y a)](l — i'(k)/e(k)) - it vanishes along 
the BZ diagonal and just off the FS when i'(k) = e(k). 
This is seen in Figs. 3, 4, and 5, which plot the ver- 




FIG. 5. The momentum dependence of the coupling con- 



stant I g(k, q 
the BZ. 



Tv/a, Tv/a) 



as a function of momentum in 



tex for q = (0, 0), (7r/a, 0), and (7r/a,7r/a), respectively. 
Thus the vertex suppresses scattering just off the FS plus 
scattering between BZ diagonals, leading to the rapid fall 
off of II with momentum transfer q as seen in Figs. (1-2). 
This is in general agreement with the results from Ref. 
Q, supporting the importance of the coupling constant. 

Finally we consider the absolute magnitude for the 
hardening of the B\ g phonon as a function of q. By 
specifying the microscopic source of the e-ph coupling 
mechanism, we can for the first time provide a numerical 
check of the strength of the coupling. For this it is useful 
to examine the neutron scattering cross section | fL5l for 
scattering by the B\ g phonon. The neutron scattering 
cross section can be written as 



dfldui 



F(q) | 2 S(q, W ), 



(12) 



where F is the form factor for scattering by the B\ g 
phonon, S = InxD (see Eq.(9)) is the phonon spectral 
function, and hu> is the difference in energy of the scat- 
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FIG. 6. The neutron cross section as a function of ui for 
different q values. The q values are chosen in equal increments 
along the diagonal (axis) of the BZ for the left (right) panel, 
respectively. Each curve has been offset for clarity. 

tered and incident neutron. As in Ref. [0, we consider 
Umklapp scattering of the neutron by the phonon for mo- 
mentum transfers of Q + q, where Q = (2n/a, 0). This 
determines the form factor as 

I F(n) I 2 - [ cos (gW 2 ) + cos(q y a/2)} 2 

1 Wl cos 2 (g x a/2) + cos 2 (q y a/2)' [ ' 

Using uJbx 9 — 348cm~ 1 and the magnitude of the Cu- 
O buckling, eE z = -0.8eV/A |§, we plot in Fig. 6 the 
results for the cross section for a d x 2_ y 2 superconductor 
for q along the diagonal and axis of the BZ, respectively. 
The e-ph coupling has been included, and the parameters 
used are the same as for Fig. (2). In addition, a constant 
phonon linewidth of 2.5 cm -1 was included which repre- 
sents the intrinsic broadening of the phonon (i.e., due to 
anharmonic lattice potentials) which is present even in 
the insulating state §. Once the spectral function was 
obtained the curve was convoluted with a Gaussian of 
width 2 meV to mimic the experimental resolution of the 
apparatus used in ||. 

Our results for the frequency shift of the phonon for dif- 
ferent momenta are compared to the experimental data 
in Fig. (7) using the above parameters. We note that the 
theory compares both qualitatively and quantitatively 
well with the data in two regards. First, the theory pre- 
dicts a hardening of the phonon for increasing momen- 
tum transfers in both directions in agreement with exper- 
iment (although data along the BZ diagonal is not very 
complete). Further, a maximum hardening of roughly 
0.9 meV observed for momentum transfer along the BZ 
axis is in agreement with the measured Q hardening of 
0.9 ± O.lmeV. The general behavior away from the zone 
center (q = 0) is reproduced as well. As mentioned the 
data along the BZ axis cannot be used to determined 
whether the energy gap changes sign, but we note that 
if the anisotropic s energy gap was used it would predict 
a larger growth of the shift along the BZ diagonals than 




FIG. 7. Comparison of the theory with the experimental 
data of Ref. §. 

displayed in Fig. (7). 

Three important details of the data are not captured 
by the theory however. The theory gives a more abrupt 
fall off for the shift along the BZ axis than seen in the ex- 
periment, although a slightly larger error bars (as seen in 
two points in Ref. would cover the discrepancy. The 
experimental data show large changes for the phonon fre- 
quency for T = 0.5T C compared to 0.1T C which the the- 
ory cannot reproduce. Lastly, if the theory is used to 
calculate the phonon frequency as a function of momen- 
tum in the normal state, the phonon develops dispersion 
that roughly tracks the dispersion seen in the supercon- 
ducting state for momentum transfers larger than na/2 
in either direction. This is also in disagreement with 
the experiment, which shows a relatively dispersionlcss 
phonon above T c . 

These points may be indicative of the neglect of strong 
electronic correlations. The theory does not include the 
effects of strong electronic correlations as the electronic 
parameters are taken in accord with local density ap- 
proximation values ||. This may crucially affect the 
momentum dependence of the lincshape as the momen- 
tum vector maps out the band dispersion as well as the 
shape of the Fermi surface. The strong local Coulomb re- 
pulsion would lead to a smaller renormalized bandwidth 
and thus a less rapid dependence of II with momentum. 
The theory could be refined once parameter choices for 
the band structure are determined via fitting to experi- 
mentally determined band dispersion and Fermi surfaces. 
This remains to be investigated. The strong tempera- 
ture dependence of the phonon even below T c remains 
puzzling. All theory curves are essential T— indepen- 
dent below T — 0.5T C as the energy gap is completely 
established. The phonon lineshape changes below T c are 
due to the opening of a phonon decay channel by break- 
ing Cooper pairs near q = and are minimized by mo- 
mentum and energy conservation at larger momentum 
transfers. Additional scattering via, e.g. impurities or 
spin fluctuations, open channels for phonon scattering at 
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small q in the normal state as well which may counteract 
the dispersion given by II. However, this cannot yield 
a strong temperature dependent scattering substantially 
below T c . 

In summary we have investigated contributions to the 
momentum dependence of the B\ g phonon resulting from 
anisotropics of the density of states, energy gap, and e- 
ph coupling constant. Using a microscopic theory for the 
origin of the e-ph coupling as given in for the first 
time the theory's prediction for the dependence of the 
lincshape of the phonon as a function of q for various 
directions in the BZ is in qualitative agreement with the 
results of Ref . Q , while the magnitude of the hardening 
of the phonon is in quantitative agreement. It is cru- 
cial that the anisotropies of all quantities are taken into 
account. 
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